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^^ ' Many of the early works in the quahty control lite rature construct cont r ol lini - 
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Wortham and Ringer! ( 19721) . 



its through the use of graphs and tables as described in 
^^ , However, the methods used in this literature are restricted to using only the values 

that the graphs and tables can provide and to the case where the parameters of the 



Jv>( ' underlying distribution are known. In this note, we briefly describe a technique which 

^ I can be used to calculate exact control limits without the use of graphs or tables. We 

also describe what are commonly referred to in the literature as fiducial limits. Fidu- 
cial limits are often used as the limits in control charting when the parameters of the 
underlying distribution are unknown. 



1 Introduction 

In many applications using attribute (count) data, the number of nonconformities often has 
discrete distributions such as the binomial, Poisson and geometric. In this note, we provide 
the fiducial limits for these distributions because, in certain cases, these limits are often 
used as the upper and lower limits for control charts. 

Early works relate d to obtaining exact control limits were based on the use of graphs and 



tables. For example, 



Wortham and Ringeii (| 19721 ) provide a good review on exact control 



limits using graphs and tables when the true parameter of the underlying distribution is 
known. These graphs and tables were used to approximate or interpolate the values needed 
to obtain an interval with 1 — a confidence level when the parameters of the underlying 
distribution were known. Before the advent of powerful computers and software, it was 
often troublesome and time-consuming to compute exact control limits. Thus, the methods 
that were based on the use of tables and graphs were prevalent in the literature because of 
their simplicity and convenience. 

Unfortunately, the studies that use these graphs and tables are restricted to using the 
specific confidence interval values that can be inferred from these tables and graphs. Addi- 
tionally, errors are difficult to avoid because of the required use of visual interpolation. How- 
ever, with powerful computers and accessible software available, it has now become quite 
trivial to calculate exact control limits without the use of graphs and tables. Exact control 
limits are obtained by calculating the corresponding quantiles of the underlying distribution 
involved. These days, there are many statistical software progr ams that provide the quan- 



tiles f or any reasonably common distribution. For example, in 



R Development Core Team 



(|2008l ). a non-commercial open source statistical software package for statistical comput- 



ing and graphics, the binomial and Poisson quantile functions are given by qbinom( ) and 
qpois( ) respectively. R can be obtained at no cost from http://www.r-project.org. 
Other commercial statistical software such as Minitab and SAS also provide the quantile 
functions for the commonly known distributions. 



Another shortcom i ng of the control hniit methods that use the graphs and tables in 
Wortham and Ringed ( 19721 ) is that, in order to use these tables, one has to assume that 
the parameter of the underlying distribution is known. In many practical examples, the 
true parameters are not known. As aforementioned, if the true parameter is known, the 
exact limits are easily calculated by using the quantile functions in the various statistical 
software packages available. If the parameters are unknown, then the control limits can 
be obtained by inverting the relation between the tails and the parameters. These control 
limits are then referred to as fiducial limits in the statistics literature and are explained in 
more detail in the following section. 



2 Fiducial Limits for the discrete distributions 



Many authors, inc 



Bickel and DoksumI (|l977l ). 



u ding 



Clopper and Pearson 



1934), [Garwood! (Il936l l. IStevensI (Il95d ) 



Kendall and StuartI (|l979l ). mention the exact or fiducial con- 



fidence limits for several discrete distributions, along with the classical confidence limits 
based on the central limit theorem. In this note, we provide a review of several fiducial 
limits which can be useful for various control charting methods including the p, u, c and g 
charts. 

2.1 The binomial distribution 



Let X denote a binomial random variable with size n and Bernoulli probability p. The 
fiducial limits for p are obtained by inverting the following two equal tails for p 



and 



EQp^(i-p)- 

A — ^. \ / 



a 
2 



a 
2' 



(1) 



(2) 



where x = 1, 2, . . . , n — 1. The lower and upper hmits with 1 — a confidence level are given 

by 
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where Foi[u^ denotes the ^ quantile of the F distribution with degrees of freedom vi and 
U2. Here, the lower limit is if x = and the upper limit is 1 if x = n. Derivations o f 
the fiducial limits for p based o n the F distribution are given by 



and 



BlvthI (119861 1. iHaldl (119521 ) 



Leemis and Trivedil (jl996l ). It should be noted that the quantile of the F distribution 



is calculated by using qf ( ) in R. 



2.2 The Poisson distribution 



Garwood! (jl936l ) provide the fiducial limits for the Poisson distribution. Let Xi,X2, . . . ,X„ 
be a random sample from a Poisson distribution with mean A. For convenience, let Y = 
Yl^=i'^j- Then, the random variable Y has a Poisson distribution with mean n\. Let 
y = y, then the lower and upper limits with 1 — a confidence level are given by 



A. = > 



(a/2) 



2n^'2 
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and 



A^ = :^r 
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where P^^ is the ^ quantile of the gamma distribution with parameters a and h. Here the 
lower limit is if y = 0. It should be noted that the quantile of the gamma distribution is 
calculated by using qgamma( ) in R. 



2.3 The geometric distribution 

Let Xi,X2, . . . ,Xn be a random sample from a geometric distribution with Bernoulli pa- 
rameter p. For convenience, let Y = X^?=i Xj. Then, the random variable Y has a negative 
binomial distribution with size n and Bernoulli parameter p. Let Y = y, then the lower 
and upper limits with 1 — a confidence level equal 
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and 
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where the upper limit is 1 if y = 0. These limits are provided in Exercise 9.22 of 



mm . 



Casella and Berger 



3 Illustrative Examples 
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